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CLOSLD  FORM  MAGNETIC  QUARTER  ORBIT  SWITCH  POINT  SOLUTION 

L.  K.  Herman 
The  Aerospace  Corporation 
£1  Segundo,  California 


Abstract 


A closed  form  solution  is  derived  for  determin- 
ing the  location  of  a set  of  magnetic  switch  points, 
spaced  90  deg  apart  in  the  orbit,  which  will  allov^ 
the  spin  axis  of  a satellite  to  be  moved  in  any, 
achievable,  desired  direction  fron^  any  initial 
orientation.  Additional  closed  forn^  equations  arc 
derived  for  determining  the  magnitude  of  the  motion 
that  will  result  for  a given  set  of  switch  points.  The 
sin^plicity  of  the  calculations  make  the  technique 
extremely  useful  for  applications  where  savings  in 
computer  time  or  storage  are  desired  and  for  design 
studies  where  quick  solutions  to  many  cases  are 
required. 

I.  Introduction 


On  a spinning  satellite,  precession  control 
torques  can  be  generated  by  the  interaction  of  a 
magnetic  coil,  which  has  its  dipole  aligned  parallel 
to  the  spin  axis,  with  the  earth's  magnetic  field. 
The  spin  axis  can  be  processed  in  any  direction  if 
the  sign  of  the  magnetic  moment  generated  by  the 
coil  is  switched  at  appropriate  points,  spaced 
90  deg  apart  in  the  orbit.  The  torque  resulting 
from  the  interaction  of  the  satellite  magnetic 
moment,  CT,.  and  the  earth's  magnetic  field,  B,  is 


(1) 


and  the  change  in  angular  momentum  resulting  from 
this  torque  is 


:: 


/ 


M X 5 dt 

s 


(2) 


This  expression  will  be  expanded  to  yield  a closed 
form  relationship  between  the  location  of  the  mag- 
netic switch  points  in  the  orbit,  the  initial  satellite 
attitude,  and  the  desired  change  in  attitude. 

11.  Magnetic  Field  Model 

The  earth's  magnetic  field  will  be  represented 
as  originating  from  a simple  dipole  rotated  11.4 
deg  from  the  geographic  pole  in  the  plane  of  the 
110  and  290  deg  east  meridians  (see  Fig.  1).  In 
the  dipole  coordinate  system,  which  has  the  z-axis 
along  the  magnetic  dipole  axis  and  the  y-axis  in  the 
plane  of  the  110  deg  meridian. 


i Mi 


0 
0 
- 1 


(3) 


where  • Ml  = 7.  95  x 10^^  pole  cm 

Transforming  to  the  nodal  coordinate  system 
(Fig.  2)  which  has  the  y-axia  normal  to  the  orbit 
plane  and  the  x-axi,  along  the  orbital  line  of  nodes 
yields. 


rigure  2.  Coordinate  frames 


(«/-90")  T^  (-6)  Ty  (4) 


where 

6 


= »gz  ♦ 


'e‘g* 


10“  - n, 

~ Greenwich  hour  angle  of  zero 


5 


earth  rotation  rate 


tg  = Greenwich  mean  time 

^n  ~ right  ascension  of  the  orbital  ascending 
node 

and  (o)  represents  the  transformation  about  the 
x«axis  through  the  rotation  angle  o,  similarly  for 
Ty  and  Tz$  i.  e.  , 


h = the  angle  between  the  x'^-axis  and  the 
x^-axis<  the  argument  of  latitude. 

P = magnetic  permeability 

The  normal  component  of  the  magnetic  field  at 
any  point  in  the  orbit  plane  is  generated  by  and 
is  oriented  along  the  -y«axis.  The  magnitude^of 
the  component  is 


T^  (o)  = 0 c(a) 


c(a)  = cos  (or) 

8(o)  s sin  (a) 

Expanding  the  equation  for  yields* 

8(11.4)  c(5) 

M"  = I M|  8f*/)s(11.4)  s(6)  + c(^)c(11.4)  (6) 

djl)  »(n.4)  8<5)  - 8(^)c(U.4) 

The  component  of  OT"  normal  to  the  orbit  plane  is 


Therefore,  the  total  magnetic  field  at  any  point 
in  the  orbit  plane,  expressed  in  the  / coordinate 
system  is 


1-2  s(e  - 0) 


B*  = K -K  /K 
n 

c (0  - e„) 


In  the  nodal  coordinate  system,  n. 


B"  = Ty  (6)  B‘ 


Expansion  yields. 


The  magnitude  of  the  component  in  the  orbit  plane 


i + (M^  )^j  (7) 

The  angle  between  the  z^-axis  and  the  component 
of  the  dipole  in  the  orbit  plane, 


\ I <.1  / 


K c(9  - 9. 


- T ■ ®m>  " T ‘<®m* 


^c(2e  - e^)  - -J-  c(0^) 


The  radial  and  ^ngential  components  of  the  magne* 
tic  field  vector  B at  any  point  in  the  orbit  plane  ^e 
generated  by  MjJ^,  In  the  / coordinate  system,  B, 

is 


-2K  s(e  - 9^) 
nn 


The  orientation  of  the  spin  axis  relative  to  the 
nodal  coordinate  system  can  be  defined  by  the  two 
Euler  angles  9 and  f (see  Fig.  3).  The  magnetic 
field  in  this  coordinate  system  is 


B®  = T^  (-7)  T^  (0)  B" 


Expansion  yields, 


B;^  c(0)  + 8(0) 

B*  = -B;^  8(0)c(T)+B^c(0)c(T)  - B^  .(T) 

-b"  8(0)  s(Y)  + b"  c(0)  s(Y)  + b"  c(T) 
* y * 


u!  M 1 
' xz- 


111.  Ma£netic  Control  Torques 

In  the  "a"  coordinate  frame,  the  magnetic  con-> 
trol  is  aligned  with  the  y-axis,  therefore, 


R s the  radial  distance  from  the  center  of 
the  earth  to  the  point  where  the  field  is 
computed. 


= M. 


®1 

0j  +2Tr 

1*  = -^  *^8  f®“<®M’^l  ‘=<2®-®m>l  (25) 

28  L -•  e. 


Substituting  the  limits  of  integration  and  switching 
the  sign  of  M,  in  each  quarter  orbit  yields 


I = -K  M 
y ns 


Evaluating  I using  the  same  technique  used  for 


lx  yields 


1 = 0 
y 


X 2 6 


^ ®(®m^  ( r -I 

Yr~  j“s[<®s+’’^2)-e3] 

+ u/2)] 

+ ^s[<®s+l''>-<®s+"'] 

-Mj(e3+2.)-(e3.|.)]j 

- '‘2®s  - ®m>] 

- M fcCZe  + 2tt  - 0 ) 

8L  8 m 

-c(2e^t.-ej] 


+ M fcize  + 3tt  - G )1 
8 1_  8 m J 


Finally* 


■'  = / 


M B " dt 
s z 


9 j + 2-it 


f ^s  / H«-®m)-<®m>]"®  I 

A ! 


Again,  evaluating  the  integral  using  the  technique 
used  for  1^  yields 


1 = M s(29  - 9 ) (33) 

z g s s m 

Substitution  of  the  integrals,  evaluated  with  Mg 
switched  at  quarter  orbit  intervals,  in  the  equation 
for  AH  ® yields, 

I s(0)  s(Y)  c(29g  - 0^)1 


6 KM  -c(Y)  s(2e„  - e„) 

s 8 m 


1 /4  orbit)  “ 9 


-c(2e  + 2s.  6^)1 

8 ^ J 


M rc(2e  + 4s  - e ) 


-c(2e  + 3s  - e ) 


I = M c(2e  - e ) (27) 

X a s s m ' ' 


c(0)  c(2eg  - ej 


If  the  sign  of  the  satellite  dipole  is  constant  over 
the  complete  orbit,  the  limits  of  integration  can, 
without  loss  of  generality,  be  replaced  with  0 and 
2s,  suid  the  integrals  become 


KM  I 

IT  jr 

■ - ®m>]| 


IT  - 6 ) 


I = / M B " dt 

y I s y 


n KM 

1 = — . s(6^) 

* 9 


-KM 

1 = S— 

y i 


-[2n  - o] 


2,  The  Batcllite  ^titude  doe«  not  change 
significantly  in  one  orbit.  Therefore, 

0 and  y are  constant 

3.  The  earth  does  not  rotate  significantly  in 
one  orbit.  Therefore, 


4.  The  spin  axis  is  constantly  aligned  with  the 
momentum  vector  due  to  damping. 

Each  of  the  integrals  will  be  evaluated  around  an 
entire  orbit,  starting  at  the  first  quarter  orbit 
switch  point,  6g,  with  Mg  positive  between  6g  and 
69  + ir/2  and  switching  signs  in  each  subsequent 
quarter  orbit  (see  Fig.  4).  The  first  term  is 
evaluated  as  follows: 


I = M 

X 8 


Figure  3.  Euler  angle  definition 

The  torque  resulting  from  the  satellite  control  coil 
interacting  with  the  earth's  magnetic  field  is, 

T*  = WT*  X (18) 


Therefore, 


de  6 dt 


-M,  ®x 


IV.  Angular  Momentum 


dt  = 4-de 
e 

Substitution  of  variables  yields 


The  change  in  angular  momentum  resulting  from 
the  magnetic  torque  is 


AH  = I T * dt 


Expansion  yields, 


AH*  = M 


-s(0)  dt  + c(0)  s(y) dt 

+c(v)  J b"  dt 

o' 

-c(«)y*B^  dt  - bw/b^  dt 


The  following  assumptions  will  be  utilized  in  evalu- 
ating the  integrals. 

1.  The  orbital  radius  and  the  orbit  rate  are 
constant.  Therefore, 

R s Constant 


■■  j’o  At  = 9 


I = M 

X 8 


j'r-" 

/ b"  ^ 

J *6 


Figure  4.  Quarter  orbit  switch  points 


-8 


Also, 


Therefore, 


KM. 


I = 


26 


j|  -ej 

- [^"  <=(  V - “]  1 


(38) 


And, 


n KM 


I = 


c (e„) 


(39) 


Substitution  of  the  integrals,  evaluated  with  M 
constant  over  the  entire  orbit,  in  the  equation  lor 
AH  * yields. 


AH 


nKM, 


(const) ' 


-8(6)  S{f)  8(0  ) 


2K 

- C(6)  sCf)-!^ 


•c(»)  c(9_) 
m 

0 


2K 


-c(0)  S(e„)  + s(0) 


(40) 


V,  Spin  Axis  Motion 


-AH“ 


AY  = 


(41) 


a6 


H c(Y) 


(42) 


where  Hq  is  the  nominal  satellite  angular 
momentum. 


VI.  Quarter  Orbit  Switch  Points 


The  quarter  orbit  switch  points,  gg,  6g  4 W2. 
and  so  forth,  can  be  determined  from  the  expres- 
sions  for  ^0  and  as  follows: 


fr  = ‘^®»  ‘ ®r 


(43) 


tan  (20 


) ( M «(tA 

^AY  c(0)  c(T)y 


(44) 


and 


^ J [^e^  + tan'*(tan(Y)sin(0)-  ^ co8(0)j 


(45) 


This  equation  defines  the  first  quarter  orbit 
switch  point  in  terms  of  the  Kuler  angles  0 and  T, 
the  angle  6n^s  which  can  be  determined  from 
Eqs.  (6),  (7),  and  (8),  and  the  ratio  of  the  desired 
change  in  attitude  A0/A'*'. 


When  the  switch  points  have  been  determined, 
the  magnitude  of  the  change  in  attitude  can  be 
determined  from  Eqs.  (41)  and  (42).  The  required 
input  data  isAh^,  Eq.  (34)  or  (40),  which  requires 
the  magnitude  of  the  satellite  dipole  and  the  orbit 
rate. 


The  units  of  the  terms  in  the  coefficients  in  the 
SfT  equations  are  discussed  in  Appendix  A.  A 
listing  of  a computer  program  for  determining 
switch  points  and  resulting  attitude  motion  is  con- 
tained in  Appendix  B. 


VII.  Performance  Evaluation 


The  spin  axis  will  be  colinear  with  the  momen- 
tum vector  if  sufficient  damping  is  present.  There- 
fore, the  change  in  the  orientation  of  the  spin  axis 
can  be  determined  if  the  change  in  angular  momen- 
tum is  known.  From  Fig.  3,  it  can  be  seen  that  for 
small  angles, 


The  performance  of  the  closed  form  equations 
was  evaluated  by  generating  switch  points  and 
attitude  change  magnitudes  for  a number  of  cases, 
using  the  computer  program  in  Appendix  B.  These 
results  were  compared  with  the  outputs  of  a sophis- 
ticated digital  simulation  using  a ninth  order  mag- 
netic field  model  while  switching  the  sign  of  the 
satellite  dipole  at  the  switch  points  determined  by 
the  closed  form  equations  and  integrating  around 
an  entire  orbit.  A summary  of  the  results  is 
shown  in  Table  1.  The  angular  differences 
between  orientation  derived  from  the  closed  form 
solution  and  the  results  of  the  more  sophisticated 
ninth  order  model  are  small  compared  to  the  level 
of  performance  generally  achieved  with  magnetic 
control  of  spin  stabilized  satellites. 


VIII.  Summary 


The  simplicity  of  the  calculations  and  the 
accuracy  of  the  results  make  this  technique 
extremely  useful  for  applications  where  savings  in 
computer  time  or  storage  are  desired  and  for 
design  studies  where  quick  solutions  to  many  cases 
may  be  required. 


Appendix  A 

Computation  of  Coefficient 


6KM 


M 

s 

L4t7  r\ 

0 

6KM 

— — = 6.  972852  x 10 

e 


,,  I M"  I M 
-23  ‘ xz'  s 


q R 


-9- 


Table  1 Performance  Evaluation 


All  angles  in  degrees  and  time  in  hours. 


1 Input  Conditions 

1 Closed  Form 

! Ninth  Order 

Angular 

Error  ♦ 

Avg 

'f  A 

Avg 

27 

SGN 

Time 

a 

s 

&0 

c 

A¥ 

c 

n 

60^Co8  t 

67^ 

d7 

0.  0 

0.  0 

1.0 

+ 

+ 

12.  0 

70.45 

1. 01 

1.01 

0.96 

1. 00 

0.  05 

0.  01 

0.  0 

0.  0 

-0.  5 

+ 

- 

18.  0 

8.41 

0.64 

-1. 28 

0.  65 

-1.23 

-0.  01 

-0.  05 

0.  0 

0.0 

2.0 

- 

- 

18.  0 

143.41 

-1.28 

-0.  64 

-1.23 

-0.62 

-0.  05 

-0.  02 

0.0 

0.  0 

0.0 

+ 

+ 

6.  5 

95.  31 

0.  00 

1. 40 

0.  03 

1.40 

-0.03 

0.  00 

0.0 

0.  0 

1.  ElO 

+ 

+ 

1.0 

43.  28 

1.  36 

0.  00 

1. 29 

-0.  05 

-0.  07 

0.  05 

30.  0 

45.  0 

-0.  33 

+ 

18.  0 

104.21 

-0.  32 

0.  97 

-0.  35 

0.  96 

0.  02 

0.  01 

30.0 

45.  0 

3.  5 

- 

- 

18.  0 

140. 91 

-1.59 

-0.  45 

-1.48 

-0.47 

-0.  07 

0.  02 

0.  0 

60.  0 

-0.  30 

+ 

- 

1.0 

6.  63 

0.  39 

- 1.  30 

0.41 

-1. 29 

-0.  01 

-0.  01 

85.  0 

0.  0 

-1.  HO 

+ 

- 

6.  5 

50.  31 

1.40 

0.  00 

1.40 

-0.  08 

0.  00 

0.  08 

0.  0 

85.  0 

-0.  3 

+ 

- 

1.  0 

6.63 

0.  39 

-1.  30 

-0.05 

-1.29 

0.  04 

-0.  01 

85.0 

85.  0 

2.0 

- 

- 

12.  0 

45.40 

-0.  02 

-0.01 

0.46 

0.  06 

0.  04 

-0.  07 

85.  0 

85.  0 

1.  ElO 

- 

- 

12.  0 

137. 95 

-1.43 

0.  00 

-1. 70 

0.  08 

0.02 

-0.  08 

♦Note:  Due  to  the  Eulerian  definition  of  0 and  f,  the  actual  angular  difference*  in  the 
0 direction,  between  the  closed  form  and  ninth  order  solutions  is  equal  to  (A0  closed 
form  - A0  ninth  order)  cos  (f  Avg).  This  effect  can  be  seen  in  the  tabular  results  when 
the  magnitude  of  the  ¥ Avg  increases. 

Orbital  Parameters  Satellite  Parameters 

R = 3761.  Onmi  Mg  = 10,000.  pole  cm 

^ = 97.71  deg  Igpin  = 183.  Slug 

0 = 0.  03  deg  ui  = 2 Rev/min 


where 


where 


6KM,  . . „ 

; ■ IS  in  ft  lb  sec 
6 


Im”  I is  in  pole  cm 
xz  ^ 


I is  in  pole  cm 
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